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We study the scaling behavior of fidehty susceptibihty density (xf) at or close to an anisotropic 
quantum critical point characterized by two different correlation length exponents u\\ and v±^ along 
parallel and perpendicular spatial directions, respectively. Our studies show that the response of 
the system due to a small change in the Hamiltonian near an anisotropic quantum critical point is 
different from that seen near an isotropic quantum critical point. In particular, for a finite system 
with linear dimension L|| (Lx) in the parallel (perpendicular) directions, the maximum value of 



with an exponent depending 
> L\'"^. We also 



for L,, 



Xi is found to increases in a power-law fashion with Z/|| for small L| 
on both and i^x and eventually crosses over to a scaling with 
propose scaling relations of heat density and defect density generated following a quench starting 
from an anisotropic quantum critical point and connect them to a generalized fidelity susceptibility. 
These predictions are verified exactly both analytically and numerically taking the example of a 
Hamiltonian showing a semi-Dirac band-crossing point. 

PACS numbers: 64.70.qj,64.70.Tg,03.75.Lm,67.85.-d 



I. INTRODUCTION 

Recent studies on fidelity and fidelity susceptibility^"— 
(xf) near a quantum critical point have contributed to a 
deeper understanding of a quantum phase transition^"— 
from the viewpoint of quantum information theory. Fi- 
delity is the measure of overlap of two neighbouring 
ground states of a quantum Hamiltonian in the parame- 
ter space. Fidelity susceptibility provides quantitatively 
the rate of change of the ground state under an in- 
finitesimal variation of the parameters of the Hamilto- 
nian. Since the ground state of a quantum many-body 
system exhibits different types of symmetries on either 
side of a quantum critical poinlj^ (QCP), a sharp drop 
of fidelity is observed right there. At the same time, 
the fidehty susceptibility usually diverges in a power law 
fashion with the system size where the exponent is given 
in terms of the quantum critical exponents^Ti^i^"— . In 
recent years a series of works have been directed to un- 
derstanding the connection between fidelity susceptibil- 
ity to quantum phase transition at critical^"— i^"— and 
multicritical points^^. Studies of fidelity per site^^^^, re- 
duced fidcliiy^'^— and geometric phase, which is also 
closely related to fidelity susceptibility^, near quantum 
critical^^i^ and multicritical^^ points have also been in- 
teresting areas of research. 

In this paper, we extend the investigation on fidelity 
susceptibility to the case of an anisotropic quantum crit- 
ical point (AQCP) and highlight the marked difference 
with the corresponding studies on an isotropic QCP. 
An interesting realization of an AQCP is seen in semi- 
Dirac band crossing point o^^'^^ where the energy gap 
scales linearly with momentum along one spatial direc- 
tion but quadratically along others unlike Dirac points 
in graphene where a gap opens linearly along both the 
directiona^i. The possibility of such a semi-Dirac point 
has been reported recentl y^^i^^ using a three unit cell 
slab of VO2 confined within insulating Ti02 and also 



in liquid He'^. A series of works on low energy prop- 
erties of a system with a semi-Dirac point has already 
been reportec&2Si^. It is to be noted that the scaling 
of defect density following a slow quench across a QCP, 
namely the Kibble Zurek scaling^"—, has also been gen- 
eralized to an AQCP using a semi-Dirac Hamiltonian^. 
An AQCP can also be realized at the edge of the gapless 
region of a two dimensional Kitaev model in a honey- 
comb lattice^i^ for which the Kibble Zurek scaling has 
also been proposed^. 

Let us consider a d-dimensional quantum mechanical 
Hamiltonian H{\) designated by a parameter A. For 
two ground state wavefunctions "00 (A) and '0o(A + (5A) in- 
finitcsimally separated in the parameter space {5\ — > 0), 
we can define fidelity {F) as^"— 



F=|(V'o(A)|Vo(A + <5A))|«l 



5X^ 



(1) 



where the fidelity susceptibility xf is the first non- 
vanishing term in the expansion of fidelity. The scaling 
behavior of xf at a QCP is well established^i^ii^. Let us 
choose the Hamiltonian to be of the form H = Hq + XHj. 
Here Hq is the Hamiltonian describing a QCP at A = 
while Hj = d\H\^_^ is the perturbation not commuting 
with Hq. One can relate the fidelity susceptibility den- 
sity (xf = l/i'^Xp) to the connected imaginary time (t) 
correlation function of the perturbation Hj(t) using the 
relation^ 



1 1 r°° 

Xt(A) = —^XF = J7il r{Hi{T)Hj{Q)),dT. 



(2) 



Using dimensional analysis in Eq. (2), we get that the 
scaling dimension of Xf is given by dim[xf] = 2A//j — 
2z + d where z is the dynamical exponent associated with 
the QCP and Ah, is the scaling dimension of the opera- 
tor Hj . Clearly a negative value of the scaling dimension 
leads to a fidelity susceptibility diverging with the sys- 
tem size L at the QCP as Xf(A = 0) - i2z-d-2AH, ^ 



2 



A positive value, on the other hand, imphes a singular 
Xf though the singular behavior appears as a subleading 
correction to a nonuniversal constanlfii. A marginal or 
relevant perturbation Hi (so that \Hi scales as the en- 
ergy) allows us to make an additional simplification com- 
ing from /S.Hj = z — \/v so that at the critical poinli^"— 



Xf 



L 



2/y-d 



(3) 



Further we get a cross-over from system size dependence 
to A dependence when the correlation length ^ ~ A^'^ 
becomes of the order of system size: 



Xf^|A|- 



d-2 



(4) 



These asymptotics are dominant for dv < 2 and sub- 
leading for di/ > 2, while at c?j/ = 2 there are additional 
logarithmic singularitie a^'^^ . 

In the following analysis, we show that the general scal- 
ing of fidelity susceptibility valid near an isotropic QCP 
gets modified due to the anisotropy in critical behavior. 
The changed scaling form naturally includes the corre- 
lation length exponents along the different spatial direc- 
tions, namely i^n and i^j^. In addition, for a finite system 
with linear dimension i|| {L±) in the parallel (perpendic- 
ular) directions, the maximum value of Xf increases with 



i|l in the limit of small L|| 



i/i^ii 
II 



< L 



) only. In 

Ly^^), wc ob- 
We also study 



contrast, for higher values of L|| (i|| 
serve a crossover and Xi scales with L_l 
the defect density, and heat density following a rapid 
quantum quench starting from an anisotropic quantum 
critical point and relate them through a generalized fi- 
delity susceptibilityii. We also highlight the connection 
to the Kibble- Zurek Scaling for the defect density follow- 
ing a slow quench through an AQCP and retrieve the 
scaling relations derived previously^i^. 

The paper is organized as follows: section II provides a 
general scaling relation of Xf associated with an AQCP. 
We do also propose the same for the heat density and 
the defect density following a general quench starting 
from the AQCP, and relate them through a generalized fi- 
delity susceptibility density. In section III we have taken 
a model Hamiltonian which shows an AQCP occurring 
in the physical systems described above and confirm our 
scaling predictions using exact analytical and numerical 
methods. Concluding remarks are presented in section 
IV. 



II. GENERAL SCALING RELATIONS 

A. Fidelity Susceptibility 

Let us consider a d dimensional quantum Hamiltonian 
showing an AQCP at A = 0. The correlations length ex- 
ponent is V — along m spatial directions and v ~ 
along rest of the {d — m) directions, called the parallel 
and perpendicular directions, respectively. The fidelity 



susceptibility as obtained from adiabatic perturbation 
theoryHiii is of the forro^ 



XF 



E 



[En - E^f 



(5) 



so that for a finite system with linear dimension L|| {L±) 
in the parallel (perpendicular) directions, the correspond- 
ing fidelity susceptibility density (xf) can be written as 



Xf 



II ^-L 



II 



{En - E^f 



Here, Eq and En denote the energy of the ground state 
and nth energy level, respectively. We can contrast 
the above equation ([6|) with the specific heat density 
(the second derivative of the ground state energy den- 
sity (£;o/i[f Li-™)) given bj^ 



XE 



1 



^d^E,/dX^ 



II ""-L 
1 



T mf d— 



En — Eq 



(7) 



Comparison of Eqs. ([6|) and ([7]) indicates that near an 
AQCP, a much stronger divergence of Xf as compared to 
Xe is expected; this is due to the higher power of energy 
difference term in the denominator of Xf ■ 

We note that near an AQCP, the specific heat xe ^ 
|A|~" where below the upper critical dimension the ex- 
ponent a satisfies a modified hyperscaling relation^ii^ii^ 
2 — a = ly^^m. + ^^{d — m) + In the limit of large 

|A| (|A| » L^'/"ii,i^'/"^), XE scales as 



XE 



(8) 



Now. in the same limit the scaling of the fidelity suscep- 
tibility density is given by^ 



1 



Xf = 



l(V'«l|fl^o)| 



XE 



^11^1"'" ^0 (^"-^o)' \En-E^\ 



_ |_^|,.||m+^_L(d-m)-2 (|A| > L||^/''ll,L-l/''-L). (9) 

In deriving Eq. (|9]) wc have used Eq. ([8]) and the fact 
that near the AQCP, En ~ Eq lAl^n^ii = jAj^^^^-^. In 
the special case oi ~ ~ v, we retrieve the expected 
scaling relation Xf ~ \'^d-2 yg^j^j^j near an isotropic quan- 
tum critical point^iii (see eq. (4)). 

On the other hand, right at the AQCP (A = 0), and 



-i/'^ii 



in the limit -^^i'"" ^ , Xf(^ = 0) scales with the 

system size L\\ as 



Xf(A = 0) 



II [lII"^' ^L^l"^). (10) 
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However, in the opposite limit L^^^'^^^ ^ ij^''^, Xii^ ~ 0) 
instead starts scaling with L±, and Eq. (fTO|) gets modified 
to 



Xf(A = 0) ~ L 



^-^"--^'-"^ (if/'^"»Li/-). (11) 



Clearly, the speeial condition L^^l'^^^ ^ lY'^^ yields 



Xf(A==0)^L 



-= — (d—m) — rn 

'"II '"II 

II 



' -jjLLm— (d—m) 



(12) 



The above scalings in Eqs. ([T0]-[T2]) suggest Xf(A = 0) 
initially increases with L|| until -^[|^''" ~ -LY'^^. Beyond 
which xt(A = 0) becomes independent of L|| and sat- 
urates to a constant value. However, in this limit the 
fidelity susceptibility density scales with Lj_ , as shown in 
Eq. dfl]). 

An alternative way of arriving at the above scalings is 
by the use of correlation functions^"— 



1 



Xf 



m r a.-ni 



LV?L 



II ^-L 



(13) 



where we define 



-Ht 



and 



with r being the imaginary time. For a relevant per- 
turbation XHj should scale as the energy, so that Hj ^ 
A''ii^i|-i. Using the relation r - L|^" and ^ A^'^n.i 

we get the scaling of Xf from Eq. (|13p given by 



^ |^|i^||m+i'i(d-m)-2^ 

which is identical to Eq. (|9]). 



(14) 



B. Heat and defect density following a sudden 
quench 

In this section we study a sudden quencb^i^ of a quan- 
tum system of amplitude A, starting from the AQCP. 
The quantities of interest are defect density^"— (riex) 
and heat density^i (Q) generated in the process. Advan- 
tage of using heat density, or the excess energy above the 
new ground state, is that it can be defined even for non- 
integrable systems. On the other hand, for an integrable 
system with non-interacting quasi-particles, it is useful to 
define defect density, which is a measure of the density 
of excited quasi-particles generated in the system. 

As A is suddenly increased from A = to its final value 
A, all the momentum modes fc|| ^ A^n and k± < A''^ get 
excited with excitation energy X'^w^w = X'^^^^ for each 
mode. This gives an excitation energy density or heat 
density of the form 



(15) 



Defect density is related to the probability of excita- 
tion, which in turn can be expressed in terms of fidelity 
susceptibilitySiiiii^. Following the above argument one 

finds that 



(16) 



Eq. (fT6| can also be derived by noticing that for a sud- 
den quench of amplitude A, all the momentum modes 
ki < A'^i I and fc2 ^ A''^ get excited with unit probability. 



givmg rie 



I rn+i'± (d—m) 



C. Generalized fidelity susceptibility density 

In this section we deal with a generic quench from an 
AQCP at time t = given by 



X{t) 



s-eit), 



(17) 



where (5 is a small parameter, and & is the step 
functionii. The case r = denotes a rapid quench of 
amplitude 5; the case r = 1 implies a slow linear quench 
with a rate 6 and so on. In all these cases the limit 5 — > 
is considered to signify a slow adiabatic time evolution. 
If the system is initially in the ground state, the tran- 
sition probability to the instantaneous excited state as 
obtained from the adiabatic perturbation theory is given 
by 



P, 



l(^»l|fl^o)P 



= <5^if;'Xl-™X2r+2, (18) 

which leads to a density of defect of the form 
1 



^11 



P 

d—m 



+2- 



(19) 



In the above, we have used the definition of a generalized 
fidelity susceptibility density xi given byii 



1 



1 



(20) 



From Eq. ([20| . one finds that xi stands for the specific 
density xe while X2 is the fidelity susceptibility density 
Xf; X4i the other hand, yields the excitation prob- 
ability following a slow linear quench starting from an 
AQCP. 

In the same spirit as in Eq. (2), a general xi can also 
be expressed in terms of time dependent connected cor- 
relation functions given by 



XI 



LJ»L'[-"(?-l) 



{Hi{r)Hi{0))cdT. (21) 



^i/^ii 



and t - - L"^ m 



Now, using A ^ i|| - ^^^^ ^ ■ - 

Eq. p7|) leads to the scaling relations L|| ^ (5 



' L 
I 



4 



^ . These suggest that one 



10 r 



5 '+""11 "I I , and T 



can further conclude Hi ~ \^\\^\\' 

(5 '+''ir'ir. Substituting for L||, i// and r 
in Eq. (1^ with / = 2?' + 2 one gets in the hmit b » 



L 



II 



X2r+2 



.i + i/^(d-ni)-2-2.^ 

i+-ii-ir 



(22) 



Therefore for a generic quench from an AQCP scahng of 
defect density gets modified to 



Up 



r+l 



6:$>L 



II 



while the corresponding heat density scales as 



.(23) 



(24) 



The expressions for Uex and Q match exactly with the 
same for fast quench (Eqs. (|15ll6p ) if we put r = 
and (5 = A, whereas, the case r = 1 correctly repro- 
duces the values for a slow linear quench starting from 
the AQCPiM3. 

The scaling relations presented above are valid as long 
as the corresponding exponents do not exceed two. Oth- 
erwise contributions from short wavelength modes be- 
come dominant and hence the low energy singularities 
associated with the critical point become subleading^iii. 



III. MODEL AND HAMILTONIAN 
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FIG. 2: Variation of Xf(A ~ 0) witli L|| as obtained numeri- 
cally for L± = 100000, 1^,1 = 1/2, i^x = 1, d = 2 and m = 1. 
Xi diverges as Xf ~ -f'li, in agreement with the scaling given 
in Eq. 
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FIG. 3: Variation of Xf(A = 0) with L± as obtained numeri- 
cally for L|| = 10000, i^ii = 1/2, = 1, rf = 2 and m = 1. Xf 
diverges as Xf ~ as expected from the scaling Eq. (|lip . 
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FIG. 1: Variation of Xf with A, as obtained numerically for 
L|l = 10000, L± = 1000, = 1/2, z/x = 1, d = 2 and m = 1. 
Xf peaks at the AQCP, and falls as |Ap'^^, as predicted in 
Eq. ©. 



We illustrate the above analytical predictions using the 
representative case of a semi-Dirac point in spatial di- 
mension d = 2. In this case, i/^ = 1/2, i/± = 1 and 



d = 2,TO = 1. In the momentum k space the Hamilto- 
nian near a semi-Dirac point can be written as the direct 
product of 2 X 2 Hamiltonians given ]3 y^6,4o,43 



A 

k 1 1 X k I 



-A 



(25) 



The fidelity susceptibility density near the semi-Dirac 
point (A = 0) can be written as as 



1 



A^l I -\- kj^ 



rdkwdk I 



./L,J.,L. (A^ + fcfl+fci)^^'^!!"'^^- 



(26) 



Rescaling fc||/-\/A = xi, kj_/ X = X2 and taking the limit 

i T 

II 



A > L|| '^,Lj^ we get 



Xf 



lAr^TT^ 7vx/AL|| A/ALx il + xi + xl) 



;dXldX2 



1 



00 poo 



dxidx2 



(27) 
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FIG. 4: Variation of Xf(-^ ~ 0) with as obtained numeri- 
cally for !/|| = 1/2, = 1, d = 2, m = 1 and L± = 10000. 
Xi saturates at > L±, as expected from Eqs. (|29l I30|l . 
Inset shows Variation of xt with Lx when L|| kept fixed at 
= 100. Xf saturates at L± > Lm. 
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FIG. 5: Kink density a function of A as obtained nu- 

merically for L|| = L± = 1000, !^|| = 1/2, i'± — 1, d — 2, 

m = 1 and A S> ii '' " , L^^''"^. riex varies as riex ~ A'^''^, as 
predicted in Eq. (|16[). 



Eq. ([29)) depends on one of the length scales, and we get 
Xf(A = 0) ~ Lu 



L 



1/2 



(for Lfi < 



(for in > i_L). 



(30) 



Numerical verifications for the scalings of Xf with L|| and 
L± discussed in Eq. (PH)) are provided in figures ([2] - U) . 
Extending our analysis of Xf to find the defect density 
following a fast quench starting from the AQCP (A = 0), 
we arrive at the scaling ~ A^Xf ~ A"^/^. This rela- 
tion is in perfect agreement with Eq. (jl6p and is verified 
numerically as shown in Fig. ([5|). However, scaling anal- 
ysis of heat density Eq. predicts Q ~ |Ap'^, which 
is subleading to the quadratic form Q ^ arising from 
contributions of short wavelength modes. This leads to 
the scaling relation Q ^ X^. This quadratic scaling is 
also checked numerically in Fig. 
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FIG. 6: Heat density Q as a function of A as obtained numer- 
ically for L|| = L± = 1000, !/|| = 1/2, z^x = 1, d = 2, m = 1 

and A ^ Lj ^ ",1/]^^''''^. Q follows the perturbative scaling 
law Q ~ A'^, as discussed in the text. 



which shows that divergence of Xf at A — >■ (see Fig. ([T])) 
and exponent (1/2), and arc in complete agreement with 
Eq. for = 2, m = 1, =1/2 and iy± = 1. 
Right at the AQCP (A = 0), we have. 



Xf(A = 0) 



1 



dfc||dfc_L. 



(28) 



The scalings k± = p/cj^ or fc| | = q^/k± simplify the above 
integral to 



1 f°° dk\\ 



dp 



tt/L± J-!T/Lii-/kT 1^ + 



(29) 



In the limit <^ L± (or ^ L±) we can approximate 
Tr/L±k?, (or 7r/L||-\/fcj_) to zero, so that the scalings in 



IV. CONCLUSIONS 

We have studied the scaling behavior of fidelity sus- 
ceptibility near an anisotropic quantum critical point. 
Anisotropic critical behaviour modifies the general scal- 
ing form of Xf- In particular, both v^^ and i'± appear in 
the scaling. In addition, even though the maximum value 
of Xf scales with L|| in the limit of small L||, at higher 
values of the same a cross-over is observed and Xf starts 
scaling with L± instead. We also propose the scaling re- 
lations for the defect density and heat density following 
a generic quantum quench starting from an AQCP and 
relate them through a generalized fidelity susceptibility. 
We have verified our general scaling predictions both nu- 
merically and analytically using the illustrative example 
of a Hamiltonian showing a semi-Dirac point. Interest- 
ingly, we show that the heat density following a rapid 
quench starting from a two-dimensional semi-Dirac point 
varies quadratically with the amplitude and the scaling 



6 



arising due to low-energy critical modes appear only as 
a sub-leading correction. 



in related works. A.D. acknowledges CSIR, New Delhi, 
for partial financial support. 



Acknowledgments 

The authors acknowledge U. Divakaran, A. 
Polkovnikov and R. R. P. Singh for collaboration 



Electronic address: lvictor@iit k.ac.inl 



Electronic address: dutta@iitk. acini 

P. Zanardi, N. Paunkovic, Phys. Rev. E, 74, 031123 (2006). 
L. C. Venuti and P. Zanardi, Phys. Rev. Lett. 99, 095701 

(2007) . 

P. Zanardi, P. Giorda and M. Cozzini, Phys. Rev. Lett. 99, 
100603 (2007). 

S.-J. Gu, Int. J. Mod. Phys. B 24, 4371 (2010). 
S. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, Cambridge, 1999). 

B. K. Chakrabarti, A. Dutta and P. Sen, Quantum Ising 
Phases and transitions in transverse Ising Models, m41 
(Springer, Heidelberg, 1996). 

M. A. Continentino, Quantum Scaling in Many-Body Sys- 
tems (World Scientific, 2001). 

V. Gritsev and A. Polkovnikov larXiv:0910.3692l published 
in "Understanding Quantum Phase Transitions", edited by 
Lincoln D. Carr (Taylor and Francis, Boca Raton, 2010). 
A. F. Albuquerque, F. Alet, C. Sire and S. Capponi, Phys. 
Rev. B 81, 064418 (2010). 

D. Schwandt, F. Alet and S. Capponi, Phys. Rev. lett., 
103, 170501 (2009). 

C. De Grandi, V. Gritsev and A. Polkovnikov, Phys. Rev. 
B 81, 012303 (2010). 

S.-J. Gu and H.-Q. Lin, Euro. Phys. Lett. 87 (2009) 10003; 
W.-L. You, Y.-W. Li, and S.-J. Gu, Phys. Rev. E 76, 
022101 (2007); S. Chen, L. Wang, S.-J. Gu, and Y. Wang, 
Phys. Rev. E 76, 061108 (2007); S.-J. Gu, H.-M. Kwok, 
W.-Q. Ning, and H.-Q. Lin, Phys. Rev. B 77, 245109 

(2008) ; S. Chen, L. Wang, S.-J. Gu, and Y. Wang, Phys. 
Rev. A 77, 032111 (2008); S. Yang, S.-J. Gu, C.-P. Sun, 
and H.-Q. Lin, Phys. Rev. A 78, 012304 (2008); W.-C. 
Yu, H.-M. Kwok, J. Gaol, and S.-J. Gu, Phys. Rev. E 80, 
021108 (2009); Z. Wang, T. Ma, S.-J. Gu, and H.-Q. Lin, 
Phys. Rev. A 81, 062350 (2010); . 

N. T. Jacobson, S. Garnerone, S. Haas, and P. Zanardi, 
Phys. Rev. B 79, 184427 (2009); A. T. Rezakhani, D. F. 
Abasto, D. A. Lidar, P. Zanardi,Phys. Rev. A 82, 012321 
(2010); S. Garnerone, D. Abasto, S. Haas, and P. Zanardi, 
Phys. Rev. A 79, 032302 (2009); S. Garnerone, D. Abasto, 
S. Haas, and P. Zanardi, Phys. Rev. Lett. 102, 057205 

(2009) ; X.-M. Lu, Z. Sun, X. Wang, and P. Zanardi, Phys. 
Rev. A 78, 032309 (2008); L. Venuti, H. Saleur, and P. 
Zanardi, Phys. Rev. B 79, 092405 (2009); D. Abasto, A. 
Hamma, and P. Zanardi, Phys. Rev. A 78, 010301 (R) 
(2008); M. Cozzini, R. lonicioiu, and P. Zanardi, Phys. 
Rev. B 76, 104420 (2007); M. Cozzini, P. Giorda, and 
P. Zanardi, Phys. Rev. B 75, 014439 (2007); L. Campos 
Venuti, M. Cozzini, P. Buonsante, F. Massel, N. Bray-Ali, 
and P. Zanardi, Phys. Rev. B 78, 115410 (2008). 

M. M. Rams and B. Damski, Phys. Rev. Lett. 106, 055701 



(2011). 

H.-Q. Zhou, and J. P. Barjaktarevic, J. Phys. A: Math. 
Theor. 41, 412001 (2008); H.-Q. Zhou, R. Ors, and G. 
Vidal, Phys. Rev. Lett. 100, 080601 (2008); B. Li, S.-H. 
Li, and H.-Q. Zhou, Phys. Rev. E 79, 060101(R) (2009); J.- 
H. Zhao and H.-Q. Zhou, Phys. Rev. B 80, 014403 (2009); 
H.-Q. Zhou, J.-H. Zhao, B. Li, J. Phys. A: Math. Theor. 
41, 492002 (2008). 
^'^ V. Mukherjee, A. Polkovnikov, and A. Dutta, Phys. Rev. 
B 83, 075118 (2011). 

H. -Q. Zhou, J. H. Zhao, and B. Li, J. Phys. A: Math. 
Theor. 41, 492002(2008); J.-H. Zhao and H.-Q. Zhou, 
Phys. Rev. B 80 014403 (2009). 

M. Znidaric and T. Prosen, J. Phys. A: Math. Gen. 36, 
2463(2003). 

2° J. Ma, L. Xu, H.-N. Xiong, and X. Wang, Phys. Rev. E 78 
051126 (2008). 

E. Eriksson and H. Johannesson, Phys. Rev. A 79 
060301 (R) (2009). 
^2 A. C. M. CaroUo and J. K. Pachos, Phys. Rev. Lett. 95, 
157203 (2005). 

S.-L. Zhu, Phys. Rev. Lett. 96 077206 (2006). 

A. Patra, V. Mukherjee and A. Dutta, J. Stat. Mech. 

(2011) P03026. 

V. Pardo and W. E. Pickett, Phys. Rev. Lett. 102, 166803 

(2009) . 

'^^ S. Banerjee, R. R. P. Singh, V. Pardo, and W. E. Pickett, 
Phys. Rev. Lett. 103, 016402 (2009). 

A. H. Castro Neto et. al. Rev. Mod. Phys. 81, 109 (2009). 

G. I. Volovik, JETP Lett. 73 (2001) 162. 

V. Pardo, and W. E. Pickett, Phys. Rev. B 81, 035111 

(2010) . 

^° P. Delplace and G. Montambaux, Phys. Rev. B 82, 035438 
(2010). 

^'^ W. H. Zurek, U. Dorner and P. ZoUer, Phys. Rev. Lett. 
95, 105701 (2005) 

A. Polkovnikov, Phys. Rev. B 72, 161201 (R) (2005). 

J. Dziarmaga, Phys. Rev. Lett. 95, 245701 (2005). 

R. W. Cherng and L. S. Levitov, Phys. Rev. A 73, 043614 

(2006). 

V. Mukherjee, U. Divakaran, A. Dutta and D. Sen, Phys. 
Rev. B 76, 174303 (2007). 

D. Sen, K. Sengupta, and S. Mondal, Phys. Rev. Lett. 101, 
016806 (2008). 

J. Dziarmaga, Advances in Physics 59, 1063 (2010). 

A. Dutta, U. Divakaran, D. Sen, B. K. Cha krabarti, T. F. 

Rosenbaum and G. Aepph. [arXiv: 1012.065 3 (2010). 

A. Polkovnikov, K. Sengupta, A. Silva, M. Vengalattore, 

arXiv: 1007.5331 (2010). 

A. Dutta, R.R.P. Singh, and U. Divakaran, EPL, 89 (2010) 
67001. 



7 



A. Kitaev, Ann. Phys. (N.Y.) 303, 2 (2003). 

K. Sengupta, D. Sen, and S. Mondal, Phys. Rev. Lett. 100, 

077204 (2008). 

T. Hikichi, S. Suzuki, and K. Sengupta, Phys. Rev. B 82, 
174305 (2010). 

G. Rigohn, G. Ortiz, and V. H. Ponce, Phys. Rev. A 78, 
052508 (2008). 

K. Binder and J.-S. Wang, J. Stat. Phys 55, 87 (1989). 



A. Dutta, B. K. Chakrabarti, and J. K. Bhattacharjee, 
Phys. Rev. B 55, 5619 (1997) 

D. Rossini, A. Silva, G. Mussardo and G. E. Santoro, Phys. 
Rev. Lett. 102, 127204 (2009). 

E. Canovi, D. Rossini, R. Fazio, G. E. Santoro and A. 
Silva, Phys. Rev. B 83, 094431 (2011). 

S.-J. Gu, Phys. Rev. E 79, 061125 (2009). 



